Many-body excitations in tunneling current spectra of a few-electron quantum dot 
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Inherent asymmetry in the tunnehng barriers of few-electron quantum dots induces intrinsically 
different tunneling currents for forward and reverse source-drain biases in the non-linear transport 
regime. Here we show that in addition to spin selection rules, overlap matrix elements between many- 
body states are crucial for the correct description of tunneling transmission through quantum dots at 
large magnetic fields. Signatures of excited {N — l)-electron states in the transport process through 
the Ai'-electron system are clearly identified in the measured transconductances. Our analysis clearly 
confirms the validity of single-electron quantum transport theory in quantum dots. 

PACS numbers: 73.21. La, 73.21.-b 



Lithographically defined quantum dots (QDs) are 
nanoscale systems in which electrons are constrained by 
quantum confinement [l|, 0] created by hetero-structure 
barriers and a spatially separated dopant charge distribu- 
tion around the dot (in the leads) modulated by external 
gate voltages. The resulting potential in the plane of 
the dot has an approximately parabolic dependence with 
distance from the center of the QD . 

Usually the QDs electronic structure is probed by 
single-electron transport spectroscopy where electrons 
tunnel through the double barrier structure under the 
influence of a source-drain Vsd 0, [1] • this approach, 
the non-linear regime, for which Vsd is of the order of 
the average energy separation between many-body states, 
gives access to the QD excited states. However, this tech- 
nique has its shortcomings since the numerous features 
observed in the current spectra Q cannot be solely ex- 
plained in terms of the QD energy spectrum. 

Indeed, during single-electron tunneling processes, the 
QD electronic configuration fiuctuates between N- and 
[N — l)-electron states, so that carrier transmission de- 
pends on the detailed occupation of the many-body en- 
ergy states [1, [3|- Hence, in a structure with asymmetric 
barriers, electrons injected in the QD through the thick 
(forward bias) rather than the thin (reverse bias) emitter 
barrier achieve different electron configurations, as in the 
former case the dot will be in the (iV — l)-electron state 
most of the time (because of the thin collector barrier 
allowing easy escape of electrons from the QD), while in 
the latter case it will be predominantly occupied by N 
electrons 

In this work we show that in a few-electron QD 
with asymmetric hetero-barriers distinct tunneling cur- 
rent spectra emerge for forward and reverse Vsd as a 
result of different non-equilibrium QD state occupations 



SHIlBl- Moreover, our analysis shows that the consid- 
eration of quantum-mechanical overlaps between many- 
body states [lol 11 1 and spin selectivity of the tunneling 
are essential for understanding the non-linear transport 
characteristics in QDs. 

Our experimental QD structure has a layered profile 
in the vertical direction, and is similar to structures used 
earlier [J, 01, «• e., the two Alo.22Gao.78As barriers and 
Ino.05Gao.95As quantum well are nominally 7.5 (thinner 
barrier with calculated tunneling rate 7^ ^ 3.6 ns~^), 9.0 
(thicker barrier with calculated tunneling rate ^ 0.5 
ns^^) and 12 nm thick, respectively. The QD formed 
in the quantum well region is controlled by a voltage 
Vc applied to a single gate wrapped around the base of 
the cylindrical mesa with a diameter of 0.54 //m [l|. The 
current- voltage characteristics are measured at a temper- 
ature ^ 0.1 K. Since the dot-lead coupling is relatively 
weak, the maximum current measured at 1.5 mV is ^ 50 
pA (conductance dl/dVsD is less than 10~^e^/h). 

Contour plots of the measured differential transcon- 
ductance dl/dVc as a function of the gate voltage Vq 
and magnetic field are shown in Fig. [Jl top row, for for- 
ward \VsD = 1-5 mV, Fig. [Ija)] and reverse [Vsd = —1.5 
mV, Fig. IHb)] source-drain biases where we take / as 
the absolute value of the current. The finite transport 
window gives rise to non-zero current over a finite range 
of Vg thus forming current stripes. Inside each stripe, 
there are somewhat weak random-looking features that 
are attributed to fluctuations in the density of states in 
the leads Q as well as a number of more pronounced 
features due to excitations in the many-electron system. 

From Fig. [U top row, one can see that for — 2, 
aside from the second stripe's lower and upper edges, 
there is only one current step within the 1.5 mV transport 
window (feature I) at low magnetic field, merging with 
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the lower edge of the stripe at about 5 T (■), followed at 
higher field by a current step (feature II) that disappears 
at - 7 T (region III). 

In order to gain insight in the current spectra, we com- 
pute the energy spectrum of the A^-electron QD by nu- 
merically diagonalizing the Hamiltonian of the system 
H M: 
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where the single-particle Hamiltonians are {h — \) 
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and the second term in Eq. ([T]) describes the Coulomb 
interaction among the electrons, m* — O.OGlrno and 
e = 13.0 are respectively the conduction band effec- 
tive mass and dielectric constant in Ino.05Gao.95As, and 
A = (i3/2)(— y, X, 0) is the vector potential in the sym- 
metric guage for a magnetic field B oriented perpendic- 
ular to the QD plane. The last term in hi accounts for 
Zeeman splitting for which we set g = —0.44. Due to 
lithographic and natural imperfections [l[, the confine- 
ment potential Vc(r) of the studied QD is known to be 
slightly elliptic, i.e., Vc(r) = (l/2)m* [ujIx^ + ^Iv^] with 
the confinement energies in the x- and y-directions being 
huj^ = 5.3 and hwy = 5.65 meV, respectively. Fig. [2] 
shows the QD energy spectra for = 2, 3 electrons. 

By comparing with the experimental spectrum in 
Fig. [2 top row, one can see that features I and II corre- 
spond respectively to the excited triplet and singlet states 
that cross at ~ 4 T 0, UBl ■ Since with increasing number 
of electrons the separation between energy levels gen- 
erally decreases, the energy window for A'^ = 3 in Fig. 
[5] contain more features whose evolution with magnetic 
field is similar to features in the N = 2 stripe. Most no- 
tably the current step (feature IV) in Fig. [T] visible only 
up to ^ 7 T (region V) is due to the lowest excited 
doublet. Strong features in the N = 3 stripe between 4 - 
6 T (Fig. [T|) can be assigned to the lowest quartet state 
crossing with the two lowest doublet states [▲ and T] 
which are clearly distinguishable despite their proximity. 
The current step at feature VII between □ (the magnetic 
field at which the singlet-triplet transition ■ occurs for 
N = 2) and ▼ due to the ground state = 3 quartet and 
the ground state N = 2 triplet 3, 17 1 is visible in both 
Fig. [T{a) and (b). The first excited spin doublet state 
is also clearly seen in Fig. [T](b) below 1 T (feature VI). 
Note that the splitting between this excited state and the 
N = 3 ground state at T due to the slight asymmetry 
of the QD is used to obtain appropriate values of lu^ and 
LUy in Vc{r). 

However, despite these agreements, there are other fea- 
tures in the experimental data of Fig. [T] which cannot be 



account for by the energy spectrum alone. For example, 
the current steps due to the lowest excited N = 2 sin- 
glet and the lowest excited = 3 doublet "disappear" 
altogether above ~ 7 T (beyond regions III and V) , even 
though the spectra calculations show these energy states 
to be present. Also, the distinctive current steps between 
O and □ (feature VII) and • and O (feature VIII) in 
Fig. [1] do not have any obvious counterparts in the cal- 
culated energy spectrum of Fig. [2] More importantly 
we can clearly see the effect of the inherent asymmetry 
of the double-barrier hetero-structure on the measured 
dl/dVc as different steps in the current at features I, II 
and VI in Figs. [Ija) and (b). 

With the eigenspectrum of the A^-electron system [Eqs. 
(U E])], we perform explicit calculations of the (sequen- 
tial) current / through the QD [lo| : 



/ = -67 ^r„^ [P„(A) + Pp{N- 1)] [Jl - fE 



(3) 
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where 7 = ^l1rI{ir + 7l) ~ 0.45 ns is the effec- 
tive "bare" electron coupling between the QD and the 
leads [ll], fL{R) = fL(R) {^Eaf} - IJ-ur)] are the Fermi 
functions that determine the energy level occupation for 
electrons tunneling in to (out of) the QD from (to) the 
left (right) lead L (R) with chemical potential hl (ma); 
and AEa/3 — Ea{N)~Ep{N ~1) is the difference between 
the a-th A^-electron and /3-th (A^ — l)-electron energies 
(the electro-chemical potential). Tap is the overlap ma- 
trix element given by 



(4) 



where \E'q(A^) and a\ are the wave function of the a-th 
A^-electron state and the operator creating an electron 
in the QD in the i-th single-particle state, respectively. 
The non-equilibrium occupation factors Pa {N) in Eq. Q 
can be determined from the steady state solution of the 
coupled master (or rate) equations [1, 0, [13] ■ In our cal- 
culations a large number of eigenstates (up to 24 for each 
A^) is included in the solution of the master equations to 
describe degeneracy effects in /S.Eap. 

The computed differential transconductance for Vsd — 
±1.5 mV is presented in Fig. [H bottom row. The maxi- 
mum calculated current is 180 pA in good agreement 
with the measured current, given mono-layer fluctuations 
in the QD structure [13] . The plots reproduce most of ob- 
served features shown in the top panels. For each stripe 
the current starts to flow (stops) when the ground state 
electro-chemical potential of the QD becomes lower than 
the chemical potential of the source /i^ (drain ^ji), as 
expected. However, within the transport window, the 
current behavior is not trivial, and depends strongly on 
the structure asymmetry (and sign of source-drain bias) 
and the magnetic field. 
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FIG. 1: Measured (top row) and computed (bottom row) differential transconductance plots for (a) the forward (Vsd = 1-5 
mV) and (b) reverse {Vsd = —1.5 mV) source-drain bias. To obtain dl/dVc in real units, the measured and computed data 
should be multiplied by 5 pA/mV and 80 pA/meV, respectively. Symbols and features (regions) I to VIII are discussed in the 
text. 



N=3 






— »^ 


S=l/2 








N=2 




^ — ^ 





3 4 5 

Magnetic Field (T) 



FIG. 2: Energy spectrum (electro-chemical potentials) for 
= 2, 3 computed with respect to the (A — l)-electron 
ground state energy within a 1.5 meV window. Dashed lines 
show the upper boundaries of the transport windows. Zee- 
man splitting of the energy levels manifests itself at higher 
magnetic fields as multiple closely separated lines unresolved 
in the experiment. Symbols and features (regions) are the 
same as in Fig. [1] 



The non-trivial behavior of the current within the sec- 
ond stripe (N — 2) in Fig. [l]can be understood by noting 
that in our strongly asymetric structure (pilHr ^ 7) 
for Vsd > 0, the QD is occupied most of the time 
by two electrons, i.e., Pa{N = 2) w ^/M, where 
M is the number of two-electron states in the trans- 
port window. On the other hand, for negative Vsd, 
Pa{N = 2) « "Yr/{2jl) 1, independent of the num- 
ber of the states. As spin-polarized triplet states for 
which \S4N) - S^N - 1)1 > 1/2 is forbidden (the over- 
lap matrix element for the spin-polarized triplet state 
Sz{N = 2) = ±1 is either zero or sa 1 depending on the 
spin of the Af = 1 system ^^(A^ = 1) = ±1/2), the cur- 
rent through the 5=1 state is always partially blocked. 
This means that with a thick collector barrier {Vsd > 0, 
feature I), the overall current decreases compared to the 
current at the onset of the stripe [l3|- However, in the 
case of a thin collector barrier {Vsd < 0), the current 
rises when the triplet enters the transport window be- 
cause in this case the Pa{N ~ 2)'s remain unchanged. 
At high magnetic fields, beyond the singlet-triplet transi- 
tion (■), the triplet becomes the ground state at the lower 
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edge of the N = 2 stripe. At the same time, the current 
increases for both positive and negative biases when the 
excited singlet state (feature II) enters the transport win- 
dow since no current bfocking is possible for this state. 
At B ~ 6 T (region III), the crossing between the two 
singlet states with different angular momenta [l^ leads 
to a decrease of the overlap integral for the lower state 
from « 0.9 to « 0.25 resulting in a reduced current step 
for this state above 7 T, and thus explaining its "fading" 
beyond region III in Fig. [TJ Note a crossing between 
the two excited doublet states in the third stripes ^ 7 
T (region V) has a similar effect on the current in that 
region. 

Features VII and VIII in the third current stripe can be 
understood in terms of single-electron tunneling j20| as 
follows. At the lower edge of the current stripe, when an 
extra electron is added to the two-electron system in the 
ground state, the three-electron ground state is formed. 
Within the stripe, however, if the transport window is 
sufficiently large, any of the three electrons in the QD 
can tunnel out leaving the two-electron system in either 
the singlet or triplet state with a probability given by 
the value of the overlap integral Q. As another elec- 
tron can be injected into the QD in any of these two- 
electron states, this process effectively gives rise to ad- 
ditional transport channels within the bias window. An 
example is feature VII between O and □ which is at- 
tributed to electron addition to the lowest N = 2 excited 
triplet forming the iV = 3 excited quartet and along 
which, in forward bias, as for feature I in Fig. [TJa), the 
current also decreases (the quartet-singlet transition is 
forbidden by the spin blockade condition [i3|)- Feature 
VIII occurs when the electro-chemical potentials AEaf^ 
associated with the lowest iV = 3 doublet and the low- 
est N = 2 triplet drop below fin. The decrease in the 
number of channels leads to a strong reduction in current 
in Fig. [Ha) and a slight increase in the current in Fig. 
[T](b). The latter happens because of the redistribution 
of electrons among available channels when the partially 
spin blocked transport channels (due to the involvement 
of the triplet state, cf. feature I) become energetically 
inaccessible [2l|. This is, however, inverse to the feature 
I in Fig. [ija) where the appearence of the corresponding 
channels in the transport window lead to the suppression 
of the current. 

In general, as the number of AEa^ with values be- 
tween and hr can be very large, these effects can 
drastically influence the QD transport properties for any 
N. However, for small Vsd quantum mechanical overlap 
^ forbids transport via the majority of the state pairs in 
the transport window (r^^ « 0), though our simulations 
do show that with increasing Vsd the current spectra 
rapidly become more complicated. For N — 2 this is ex- 
pected to occur for biases exceeding 2 mV in the magnetic 
field range of interest when the first excited one-electron 
state pops into the transport window. Also, at high mag- 



netic fields, the large density of states and the number of 
the transport mechanisms involved will result in smaller 
changes in current for each individual transport chan- 
nel [due to the relatively small variations among Pq(A^)], 
leading to overall decrease in visibility of features in the 
current stripes. 

In the above analysis, we assumed no momentum or 
spin relaxation (2^ . Since the slowest tunneling time in 
our QD is considerably shorter than the spin relaxation 
time which has been measured to be ~ 100 /zs in verti- 
cal QDs [isl, the neglect of spin relaxation should not 
affect noticeably the current, i.e., electrons will tunnel 
out of the QD long before the relaxation causes a change 
in the level occupation. Also, assuming the momentum 
relaxation time ~ 10 ns measured for = 1 23| remains 
valid for A'^ > 1, this relaxation mechanism may lead to 
a partial suppression of the current steps for some of the 
excited states for Vsd > 0. For example, this could ex- 
plain why the first excited A^ = 3 spin doublet, already 
expected to be rather weak, is not noticeable (no fea- 
ture VI) at low magnetic fields (< 1 T) in Fig. [Ha), top 
row [13. 
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